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We explore the spontaneous formation of an excitonic insulator state at the semimetal- 
semiconductor transition of mixed-valence materials in the framework of the spinless Falicov-Kimball 
model with direct /-/ electron hopping. Adapting the projector-based renormalization method, we 
obtain a set of renormalization differential equations for the extended Falicov-Kimball model param- 
eters and finally derive analytical expressions for the order parameter, as well as for the renormal- 
ized c- and /-electron dispersions, momentum distributions, and wave-vector resolved single-particle 
spectral functions. Our numerical results proved the valence transition picture, related to the ap- 
pearance of the excitonic insulator phase, in the case of overlapping c and / bands. Thereby the pho- 
toemission spectra show significant differences between the weak-to-intermediate and intermediate- 
to-strong Coulomb attraction regimes, indicating a BCS-BEC transition of the excitonic condensate. 

PACS numbers: 71.28.-(-d, 71.35.Lk, 71.30.-fh 



I. INTRODUCTION 

The idea that an excitonic phase appears — under 
certain circumstances — at the semiconductor-semimetal 
transition dates back about half a century^'^. The forma- 
tion of excitons is driven by the Coulomb attraction be- 
tween conduction-band electrons and valence-band holes. 
Provided a large enough number of sufficiently long-lived 
excitons was created, a subsequent spontaneous conden- 
sation of these composite Bose quasiparticles may set in. 
The excitonic instability is expected to happen, when 
semimetals with very small band overlap or semiconduc- 
tors with very small band gap are cooled to extremely 
low temperatures^. The excitonic condensate typifies a 
macroscopic phase-coherent, insulating state, which sep- 
arates the semimetal from the semiconductor (see Fig.[T]). 

Surprisingly, to date, there is no free of doubt re- 
alization of the excitonic insulator (EI) state in na- 
ture. Nowadays experiments report data, however, which 
strongly support the theoretical predictions of the EI 
phase. Along this line, experiments on coupled quantum- 
well structures, e.g., have shown unusual properties 
which were inferred as indications of excitonic conden- 
sation^. Temperature dependent angle-resolved photo- 
electron spectroscopy (ARPES) on lT~-TiSe2 transition- 
metal dichalcogenides are in favor of the EI scenario as 
driving force for the observed charge-density-wave tran- 
sition^. X-ray photoemission spectroscopy and ARPES 
on quasi one-dimensional (ID) Ta2NiSe5 reveal that the 
ground state can be viewed as EI state between the 
Ni 3d-Se Ap hole and the Ta 5d electron''. Further 
real-system candidates for the EI state are pressure- 
sensitive rare-earth chalcogenides, such as mixed- valence 
TmSeo.45Teo.55. For this compounds electrical and ther- 
mal (transport) measurements indicate exciton conden- 
sation, at temperatures below 20 K in the pressure range 
between 5 and 11 kbar-. 

Also from the theoretical side the existence of the EI 



is still controversial. Most of the early mean-field ap- 
proaches work with an effective- mass Mott-Wannier-type 
exciton model and exploit the analogy to the BCS theory 
of superconductivityiS (for a more recent calculation of 
the phase diagram see Refs. lolfTll). Here the major prob- 
lem is that the excitonic phases (excitonic gas, EI) turn 
out to be unstable against a metallic electron-hole liq- 
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FIG. 1: (Color onhne) EI formation and BCS-BEC transi- 
tion scenario. At the semimetal-semiconductor transition the 
ground state of the system may become unstable with respect 
to the spontaneous formation of excitons near the point at 
which band overlap occurs. Starting from a semimetal with 
small density of electrons and holes (such that the Coulomb 
interaction is basically unscreened), the number of free carri- 
ers varies discontinuously under an applied perturbation, sig- 
naling a phase transition^ . Approaching the transition from 
the semiconductor side, an anomaly occurs when the band 
gap, tuned, e.g., by external pressure, becomes less than the 
exciton binding energy^. Depending on from which side of 
the semimetal-semiconductor transition the EI is reached, the 
EI can be viewed either as BCS condensate of loosely-bound 
electron-hole pairs or as BEC of preformed tightly-bound exci- 
tons- . A finite order parameter A indicates the new distorted 
phase of the crystal, with coherence between conduction- and 
valence-band electrons and a gap for charge excitations. 
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yj(|i2,i3^ At present, Falicov-Kimball-type models seem 
to be the most promising candidates for realizing collec- 
tive exciton phases. This particularly holds for the ex- 
tended Falicov-Kimball model (EFKM), which includes 
a direct /-/ electron hopping term, that — having again 
the Tm[Se,Te] system in mind — is certainly more real- 
istic than entirely localized / electrons. By means of 
unbiased constrained path Monte Carlo (CPMC) sim- 
ulations the EFKM has been proven to exhibit criti- 
cal excitonic correlations (an EI ground state) in case 
of ID (2D}iiii^. Subsequent Hartree-Fock calculations 
yield the ground-state phase diagram of the 2D EFKM 
in excellent agreement with the CPMC data^, support- 
ing the applicability of such mean-field approaches also 
in the 3D situation^'^'i'^. For the 3D EFKM, the existence 
of the EI phase was corroborated by more sophisticated 
slave-boson approaches^^iiS. 

Assuming that f-c electron coherence may lead to an 
EI phase in the EFKM, the properties of the excitonic 
state should be explored in more detail. In this regard, 
the anticipated 'BCS-BEC crossover' scenario^iiii^, con- 
necting the physics of BCS superconductivity with that 
of Bose-Einstein condensates (BECs), is of vital impor- 
tance. Calculating the Frenkel-type exciton propagator 
within a random phase approximation scheme, the exis- 
tence of excitonic bound-states has been established for 
the EFKM also above Tc, on the semiconductor side of 
the semiconductor-semimetal transition^. No bound- 
states were found on the semi-metallic side. Accord- 
ingly the condensation process should differ by its na- 
ture: While 'exciton' formation and condensation simul- 
taneously take place on the semi-metallic (BCS) side, 
preformed excitons will condense on the semiconducting 
(BEC) side as the temperature is lowered (cf. Fig. [1]). 
Looking at the EI order parameter A only, it seems diffi- 
cult to examine the BCS-BEC crossover. The gap equa- 
tion for A will of course not discriminate between both 
regimes. Photoemission spectroscopy, on the other hand, 
will directly probe the elementary excitations and energy 
dispersion and therefore provides extremely useful infor- 
mation about the BCS-BEC crossover. This has been 
shown quite recently in the context of ultracold (atomic) 
Fermi gases^. 

In this work, we will follow this perspective and ex- 
amine the EI phase in terms of the EFKM particularly 
with regard to a BCS-BEC crossover. To this end, we 
analyze the equilibrium and spectral properties of the 
model at zero temperature, using the so-called projective 
renormalization method (PRM)^^'^'*. This technique has 
already been successively applied to a great variety of 
many-body problem o^^'^^ . Here we calculate and discuss 
the photoemission spectra of the EFKM in order to probe 
the signatures of the excitonic condensate. The paper is 
organized as follows. Section II introduces the EFKM. 
The theoretical approach is outlined in Sec. Ill, where 
the general concept of the PRM is resumed in Sec. III. A 
and explicit expressions for the renormalization differen- 
tial equations, particle number expectation values, cor- 



relation functions and single-particle spectral functions 
are given in Sec. III. B. Section IV presents the corre- 
sponding numerical results. Our main conclusions can 
be found in Sec. V. 



II. EXTENDED FALICOV-KIMBALL MODEL 

The Hamiltonian for the EFKM is written 

^ = E ^"k^Ck + E ^"k/k/k + E f^^?"/ ' (1) 

k k i 

where cj^. (c^.) and fl (f^) are the creation (annihila- 
tion) operators in momentum (k-) space of spinless c- 
and /-electrons, respectively, and and n{ are the 
corresponding occupation numbers in real space. The 
Fourier-transformed fermionic operators are defined via 
= -ij^'^j^Tije''^^^ , where rj = c,f, and the 77-fermion 
dispersion is 

el = - t^Tk - M (2) 

with on-site energy e^. In Eq. fjL denotes the chemical 
potential. In the tight-binding limit, on a D-dimensional 
hypercubic lattice, we have 7^ = 2^^^^cosfcrf. The 
sign of f^t^ determines whether we deal with a direct 
[f^t^ < 0) or indirect {f't-^ > 0) band-gap situation. 
Usually, the c-electrons are considered to be 'light' and 
their hopping integral is taken to be the unit of energy 
(t'^ = 1), while the /-electrons are 'heavy', i.e., \t^\ < 1. 
For 1^=0 (dispersionless / band), the local /-electron 
number is strictly conservedi^. The third term in the 
Hamiltonian ([1]) represents the Coulomb interaction be- 
tween c and / electrons at the same lattice site. Hence, 
if the c and / bands are degenerate, e'^ = and = , 
the EFKM reduces to the standard Hubbard model. 

In order to address the formation of the EI state in the 
EFKM, we look for a non-vanishing excitonic expecta- 
tion value (c^/), indicating a kind of spontaneous sym- 
metry breaking due the pairing of c electrons {f^ > 0) 
with / holes {t^ < 0). This is quite similar to the prob- 
lem of electronic ferroelectricity, where {c' f) ^ causes 
electrical polarizability without an interband transition 
driving field, provided the c and / states have different 
parityi^iii. Thereby, depending on the sign of i/ (direct 
or indirect gap), ferro- or antiferroelectric phases may 
exist. 

To proceed, we introduce two-particle interaction op- 
erators in momentum space, 

"■kikska = CkjCk^/k^/k^+kg-ka ' (3) 

and rewrite the EFKM Hamiltonian ([1]) in a normal- 
ordered formal: 

^ = E^k:4ck:+E4:/k/k: (4) 

k k 

- E ^ /k^k : +H.C.) + ^ E ^ «k,k.k3 : , 

k kik2k3 
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where 

k 

with dp. ~ (cj(;/k)i plays the role of the EI order param- 
eter. Note that choosing the normal-ordered represen- 
tation of operators, the symmetry of the Hamiltonian is 
explicitly broken, and iterating the self-consistency equa- 
tion derived below will readily give (meta-) stable solu- 
tions^. In the Hamiltonian the on-site energies were 
shifted by a Hartree term, 

^c(/) ^ _C(/) ^ ^^^/(,)^ ^ 

where (n'') ~ ^kiv^Vk) the particle number den- 
sities of c or / electrons for a system with lattice sites. 
In what follows, we consider the half- filled band case, i.e., 
we fix the total electron density n — {n'^) + {n^) — 1. 

III. THEORETICAL APPROACH 
A. Projector-based renormalization method 

The PRM was recently developed with the aim to di- 
agonalize many-particle systems^. One of the main ad- 
vantages of the method is to find broken symmetry so- 
lutions of phase transitions^. For example, in Ref. 
the method was successfully applied to the t — J model 
in order to study superconducting d-wave solutions for 
cuprates. 

The PRM starts from the decomposition of a given 
many-particle Hamiltonian into an 'unperturbed' part 
T-Iq and into a 'perturbation' Tii, where the unperturbed 
part Tio should be solvable. Suppose all diagonal matrix 
elements of Jii between eigenvectors of T-Lq vanish, the 
part Til accounts for all transitions between the eigen- 
states of Ho with nonzero transition energies. Then, the 
first goal of the PRM is to transform the initial Hamil- 
tonian into an effective Hamiltonian Hx which contains 
no longer transition operators with energies larger than 
some chosen cutoff A. Thereby, the Hamiltonian Tlx is 
formally obtained by applying a unitary transformation 

-Ux^e^'He-^\ (7) 

The transformed Hamiltonian Ha, which has the same 
eigenspectrum as the original Hamiltonian T-L, can again 
be decomposed into two parts 

nx=Ho,x+Hi,x. (8) 

Due to construction, all matrix elements (tiaIHi.aI^t-a) 
of "Hi, A with energy differences \E}l — > A should 
vanish, i.e. (nxlHi^xliTix) = 0, where and \nx) are the 
new renormalized eigenvalues and eigenstates of ?^o,a- 
Note that neither \nx) nor \mx) have to be low-energy 
eigenstates of Ho, a- To ensure hermiticity of "Ha, the 



generator Xx of the unitary transformation has to satisfy 
^1 = -'^^- . 

A crucial idea for the elimination procedure in the 
PRM is to introduce generalized projection operators Pa 
and Qa = 1 — Pa- Here Pa is defined by 

Pa^= J2 \nx){mx\{nx\A\mx), (9) 

|_E,-^-£;^J<A 

applied on any operator variable A of the Hilbert space of 
the system. Note that in expression ^ only states \nx) 
and \mx) satisfying \E^ — E^\ < A contribute to the 
transition matrix. Thus, Pa projects on the low energy 
transitions of A, whereas the orthogonal complement Qa 
projects on the high-energy transitions of A. To find an 
appropriate generator Xx for the unitary transformation 
from H to Hx, the obvious relation 

QaHa = (10) 

has to be fulfilled. 

In the original version of the PRM23, the ehmina- 
tion procedure is performed step-wise. Suppose A is the 
largest transition energy of the original Hamiltonian "H, 
in the first elimination step all transitions in an energy 
shell of width AA between A and A — A A will be removed. 
The subsequent steps remove, roughly speaking, all tran- 
sitions in the next shell of width AA between A — AA and 
A — 2AA, and so on. The unitary transformation for the 
intermediate step from a cutoff A to the new cutoff A— A A 
reads 

Ha-aa = e-^"'^^ Hx e-""''^' . (11) 

Here, the generator Xx^ax has to fulfill the requirement 

Qa-aaHa-aa - 0, (12) 

in analogy to (|10p . Thus Hx-AX has no matrix elements 
that connect eigenstates of Ho,a-aa with energy differ- 
ences larger than A — AA. 

By help of Eqs. (fTT|) and the generator Aa^aa can 
easily be constructed in perturbation theory with respect 
to Hi, A- Up to first order it reads 

--fA.AA = Qa-AaHi.a = Qa-AA-^A.AA • (13) 

Lo,A 

Here, Lq.a is the Liouville superoperator of the unper- 
turbed Hamiltonian Ho, a, which is defined by the com- 
mutator of Ho, A with any operator A on which Lo,a is 
applied, i.e. 'Lo,xA = [Ho,a,^]- Note, however, that 
the generator A'a,aa is not completely fixed by Eqs. (fTTj) 
and In fact, only the part Qa-aa-^a.aa is deter- 

mined by equation (IT^ . The part Pa-aa-^a.aa with only 
low-energy transitions can still be chosen arbitrarily and 
was set identical to zero in Eq. (jl3|) . Any physical quan- 
tities, which is evaluated in the framework of the PRM, 
is independent of a particular choice of Pa-aA'%'a,aa'^- 
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This freedom can be used to derive a continuous ver- 
sion of the method. Thereby, the low excitation part 
Pa-aa-^a.aa is chosen proportional to AA, which al- 
lows to derive differential equations for the A-dependence 
of the parameter values in the Hamiltonian during the 
rcnormalization procedure. As in the discrete versioi*^, 
also in the continuous version the elimination starts from 
the original model and proceeds until A = 0. At this 
point, all transitions operators from "Hi have been used 
up and the final Hamiltonian is diagonal or at least quasi- 
diagonal which allows to evaluate expectation values. 
Note that the parameters of the renormalized Hamil- 
tonian depend on the parameter values of the original 
model %. 

To evaluate expectation values of operators formed 
with the full Hamiltonian, we have to apply the unitary 
transformation to A as well, 

{A) = -T^-^ - {A{\))h. = {A)^ , (14) 

where we define A{X) = e'^'-Ae^'^'', A = A{X -J> 0), and 
Ti, = Ha->-o. Thus additional rcnormalization equations 
are required for ^(A). 



B. Application to the EFKM 

1. Renormalization equations 

In order to derive the renormalization equations for the 
parameters of the Hamiltonian, we first decompose the 
original Hamiltonian H into two parts 

H = no+Hi, (15) 

where 

k k 

+ E ^ : +H.C.) , (16) 

k 

and 

^1 = ^ E ^«k,k.k3: • (17) 

kikska 

Note again that the perturbation Hi only contains the 
fluctuating operator part of the Coulomb repulsion oc 
U. Following the ideas of the PRM approach, we make 
the following ansatz for the renormalized Hamiltonian 
Hx after all transitions with energies larger than A are 
integrated out: 

Ha=Ho,a+Hi,a (18) 



with 

^o.A = E=k,A : 4ck : +E=k,A : /k/k : +E>. 

k k 

+ E(Ak,A:/^k:+H.c.) , (19) 

k 

Hi,A = ^Pa E t^k,k.k3.A : «k,k.k3 : • (20) 

kik2k3 

Here, Pa projects on all low-energy transitions with re- 
spect to the unperturbed Hamiltonian Ho, a which are 
smaller than A. Due to renormalization all prefactors in 
Eqs. (fT9|). (|20l) may now depend on the wave vector k 
and on the energy cutoff A. The quantity E\ is an energy 
shift which enters during the renormalization procedure. 
In order to evaluate the action of the superoperator on 
the interaction operator in T-Li^\ one has to decompose the 
fluctuation operators : atj^k^ka ■ i^^^o eigenmodes of Ho, a- 
Obviously, the diagonalization of "Ho, a requires an addi- 
tional unitary transformation. However, for the values 
of U, used in the numerical evaluation below, the mixing 
parameter Ak,A in Eq- turns out to be always small 
compared to the energy difference kk a ~ ^k aI- "^^^^ fol- 
lows from the Hartree shifts of the one-particle energies 
in Eq. ([6|). Thus, using as approximation Lo,ac|^ = ^k'^k 
and Lo,A/k = ^k/k' '^^^ conclude 

Hi, A = E ^('^ ~ l%ik2k3,Al) C^kik2k3,A • '^kikaka ' ' 

kik2k3 

(21) 

where 

%ik2k3,A = £ki,A - Ska, A + ^LaA ~ '^ki+k3-k2,A (^2) 

is the approximate excitation energy of : akik2k3 ■' 

Lo,A : akik2k3 ■= %ik2k3,A '■ "kikaka ' ' (23) 

The 0-function in Eq. (I21[) ensures that only transitions 
with excitation energies smaller than A remain in Hi,a- 

By integrating out all transitions between the cutoff 
A of the original model and A = 0, all parameters of the 
original model will become renormalized. To find their A- 
dependence, we derive renormalization equations for the 

parameters ej^ ;)^, ejj^ ^, Ak,A, and t/kik2k3 a- The initial 
parameter values are determined by the original model 
(A = A): 

4,A = 4, Ak^A = ^, (24) 

Sk,A ^ ^k ' f^kik2k3,A = • (25) 

Note that the energy shift E\ in Ho, a has no effect on 
expectation values and will be left out in what follows. 

Next we have to construct the generator Xx.AX of 
transformation pT|) . Using relation the high transi- 
tion energy part reads in lowest order perturbation the- 
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ory according to Eqs. p3l 



Qa-aa-^a.aa 



- T 



kik2k3,A 



1-e 



kik2k3,A-AAj 



^ kik2k3 '^kikaka^A 

x8kik2k3,A : akik2k3 • ' (26) 



where we have defined 9kik2k3,A = 6(A - |j?kik2k3,Al)- 
Here the product of the two 0-function assures that only 
excitations between A — AA and A are ehminated by the 
unitary transformation (1111) . As mentioned before, in the 
present approach we prefer to use a continuous version 
of the PRM approach which is based on the choice of the 
orthogonal complement part Pa-aa-^a.aa of the gener- 
ator. Thereby, Pa-aA'^a,aa is chosen proportional to 
AA, which means that Qa-aA'^a,aa can be neglected in 
the limit AA 02i. With A'a,aa « Pa-aa-^a.aa, the 
following operator form for the generator can be used 



X: 



AA 



k2k3,A 0kik2k3,A-AA 



kik2k3 



x©kik2k3,A : akik2k3 • ' (27) 



A.AA 



where the operators are taken over from expression 
Note that the two 0-functions guarantee that expression 
(1271) corresponds to the generator part with low energy 
excitations only. For the coefficients Q:kik2k3 a ^® make 
the following ansatz: 



'^kik2k3,A 



'7klk2k3,A 



k(A 



|?7kik2k3,A|)2 



kik2k3,A 



(28) 



which is an appropriate choice in the continuous version 
of the PRIVP*. The constant n in (^5]) denotes an en- 
ergy constant to ensure that the parameter Q!kik2k3 a 
has the correct dimension of an inverse energy. At 
first glance, one might expect that Q!kik2k3,A diverges at 

= |^kik2k3.A|- Instead it vanishes exponentially at this 
point which follows from the renormalization equation 
for C^kik2k3,A' given below. 

Our aim is to derive renormalization equations for 
renormalized Hamiltonian. The transformation (|lip re- 
lates the Hamiltonian T-Lx at cutoff A to that at the re- 
duced cutoff A — AA. With Eq. (P7| one finds in the limit 
AA -> 0: 



^ X! '5kik2k3,A0kik2k3,A[: akik2k3 -'^^ 



dUx 
d\ 

(29) 

Note that the evaluation of the commutator also leads to 
new operators which are not present in the ansatz (|18p for 
T-Lx- Therefore an additional factorization has to be used. 
Comparison with the generic derivation of (fT8|) leads to 
the following set of coupled renormalization equations 
which describe the A-dependent renormalization of the 
parameters of Hx'- 



de 



k,A 



^2 -'kikk2,A"kki,ki-|-k2-k 
kik2 

X((<+k2-k)-(<)) 

$Z ^kkik2,AO;kik,k+k2-ki,A("-ki 



a(1-«)) 



kik2 



X(("(+k2-ki )-(<)) 



(30) 



de 



k,A 



dX 



X/ ^kik2,k-ki+k2,A"k2kik,A("k-ki+k2) 
kik2 

X(«)-«2)) 

X! ^kik2k,A"k2ki,k+ki- 



-k2,A 



kik2 



dA 



k,A 



d\ 



x(l-("k+k,-k2))(«) (31) 

-T^5I"kikk,AAk^ A(("ki) - ("ki)) 



ki 



{t^kik2k2,Aakkiki,A'^k(l - 



kik2 

^kikk,A'5kk2k2,A^ki (""k) 



kik2k,A 



ki+k-k2,kik2,A^ki+k-k2("k) 

+ "k2k,ki+k-k2,A4i(l - (4i+k-k2)) }■ 

(32) 

Here, we have defined expectation values 

«) = (4ck), {nl) = (/^/k), dk = (4/k), (33) 

which are formed with the full Hamiltonian. There is 
also an additional renormalization equation for the A- 
dependent coupling C/kik2kA- 1* reads 



dU, 



kik2k,A 

~dX 



'7kik2k,A"^kik2k,A • 



(34) 



Integrating the whole set of differential equations with 
the initial values given by Eq. (|25|) , the completely renor- 
malized Hamiltonian % := Hx^o ~ ^o,A-i.o is obtained 

^ = E=~k:4ck:+E4^/k/k^ 

k k 

+ ^(Ak:/K:+H.c.), (35) 

k 

where the quantities with tilde sign denote the parameter 
values at A 0. The final Hamiltonian (|35l) can be 
diagonalized by use of a Bogoliubov transformational: 

^ = E^^4^k:+E^k:/k/k:+^- (36) 
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Here, and are the new quasiparticle operators 



with 



/k = -^^k4 + "k/k: 



(37) 
(38) 



1 



1 - sgn(4 - Ek) 



-/ ~c ' 



The quasiparticle energies are given by 
£k + € sgn(e^ - e^) 



(39) 
(40) 



E, 



2 2 
/ _ £k + gk , sgn(gk-gk) 



t^k , (41) 



t^k , (42) 



where 



2. Expectation values 



(43) 



The expectation values ((33|) in the set of renormaliza- 
tion equations have to be evaluated self-consistently. Ac- 
cording to relation (1141) , thereby the same unitary trans- 
formation as for the Hamiltonian has to be used. For 
instance, following Eq. (fT4|) . the expectation value (n^) 
can be expressed by 



('^k) - (4ck) - (4(A ^ 0)ck(A 0))^ , 



(44) 



where the average on the r.h.s. is formed with the 
fully renormalized Hamiltonian H, and c|j.(A) is given by 
c|j.(A) = e'^^cjj.e"'^^. For the transformed operator we 

use as ansatz 



4(A) 



'k,A 



X! 2/kikk2,A ■ Cki/k2'^ki-(-k2-k ■ ■ (45) 
kika 



with a coherent part oc x^^ ^ and an incoherent part 
DC ykikk2 A- The operator structure in (145]) is again taken 
over from the lowest order expansion of the unitary trans- 
formation. For the A-dependent coefficients and 
2/kikk2 A renormalization equations can be derived. 
They read 



dx 



k,A 



-y 



(IX ]\[2 2/kikk2,A"kki,ki+k2-k,A 

kik2 



X (("C+k2-k) - (<)) + - (<+k2-k)) 

'^ykikk2,A _ 



dX 



''k.A'-'kikka.A- 



Integration between A (where = 1 and 2/kikk2 A ~ 0) 
and A = leads to 

: 4(A ^ 0) := £k : 4 : 

+ X/ ^kikk2 ■ '^ki'^k2-/'ki+k2-k • ' (48) 



kik2 



from which (n^) is found 

^ a72 X! I^kikk2| (CkiCki)«(/k2/k2)?i 



kik2 



^1 - (/k,+k2-k/k.+k2-k)^) • (49) 

The remaining expectation values (n^) and dj^ can be 
evaluated by using an equivalent ansatz for : /k('^) '■, 

■ flW •= ^k,A • /k • 

^ ]^ ^ 2/kik2,k-ki+k2,A • Cki'^k2/k-ki+k2 



(50) 



kik2 



consisting again of a coherent and an incoherent part 
with A-dependent coefficients x'-^ ^ and 2/kj^k2 k-ki-i-k2 Aj 
respectively. Their renormalization equations read 



dx 



k,A 



d\ 



1 



k2,k-ki+k2,AQ!k2kik,A 



("-k-ki+ka) 



kik2 



x(«)-«2)) + ("k2)(l-«))j. (51) 

'^ykik2,k-ki+k2,A _ ~ ^j.2-j 



dA 



''k,A"kik2,k-ki-(-k2,A 



where the initial values are 

2^kik2,k-ki+k2,A ^ 

led to 

{ni) = K\'{flh)n 



= 1 and 
Similar to Eq. we are 



f-k.A 



N2 E I?^kik2,k-ki+k2 I (/k-ki+ki /k-ki+k2 )« 



kik2 



4iCki)« (l - (42Ck2)w 



(53) 



^k — ^k^k(/kCk)-H 

" a72 E ykik2,k-ki+k2j'kik,k-ki+k2(CkiCki)« 



kik2 

ft 



(54) 



x(/k2'^k2)«(/k-ki+k2/k-ki+k2)w • 

In the last step, one has to evaluate the expectation val- 
ues on the r.h.sides of Eqs. (gH), ((531), and (HH), which 
are formed with H. Using the diagonal form of T-L in ([55)1 

one easily finds'^^ 

(55) 
(56) 



))", 

(46) 


('^k'^k)-fi 




^vlfiEl), 


(47) 




= -[fiK)- 


- f{Ei)]sgn{ei 



7 



where f{E^) is the Fermi function. 



C. Spectral functions 

Let us consider the one-particle spectral function for c 
electrons 



1. 



A''{k,uj) = ImG"=(k,cj), 



(58) 



where G'^(k, w) is the Fourier transform of the retarded 
Green function 



G^(k,w) = ((ck;4))(c. + *0+) 



(59) 



Using again relation (fM)) . the Green function can be 
rewritten 

((ck; cl)) (a.) = ((ck(A 0); 4(A ^ 0)))^(w) , (60) 

where the expectation value on the r.h.side is again 
formed with T-L. Using expression dUl) for cliX -■ 0), we 
are immediately led the following result for the c-electron 
spectral function 



kika 



[("^ki'^ki)?i ((/k2-/'k2)'H (•/'ki+ka-k/ki+ka-k)?^ 



(/ki+k2-k/ki+k2-k)?i (l (/k2-^k2)-H 



(61) 



Note that we have restricted ourselves to the leading or- 
der in the EI order parameter. In the same way, we can 
also calculate the spectral function A^ik^uj) for the /- 
electrons. The final result reads 



A/(k, 



|2 



J^p E l^ki 



vlS{u;-E^)+ulS{uj~El) 

|2 



k2,k-ki+k2 



kik2 



'^[^-(^k^-k.+k2-^k2+^kj] 
(cLckJ?i (l - (Ck2Ck2)«) 

+ ((42Ck2>?i " {(i,Ck,)n} 

^ ^ (/k-ki+k2/k-ki-|-k2)?i 



(62) 



IV. NUMERICAL RESULTS AND DISCUSSION 

We now will evaluate the analytical expressions of the 
PRM approach outlined so far. Of course, the set of 
equations (1551) . and ([M)) has to be solved numeri- 

cally. To this end, we choose some initial values for (n^), 



{nl), and d^. (assuming {cl^fi^) = (/kCk))' ^^'^ determine 
the renormalization of the Hamiltonian and all operators, 
by solving the differential equations (|30|) -(|32 |) . (|47|) . and 
([5^. Performing the limit A 0, all model parame- 
ters will be renormalized. Then, using the new ex- 
pectation values [Eqs. ([55l) -(|57 |) ] are calculated, and the 
renormalization process of the Hamiltonian is restarted. 
Convergence is assumed to be achieved if all quantities 
are determined with a relative error less than 10^*^. The 
dynamical correlation functions (|5T|) -- (I5^ are evaluated 
using a Gaussian broadening in energy space of width 
0.06. Because of the large number of differential equa- 
tions that have to be solved, we confine ourselves, in what 
follows, to the investigation of the ID case and limit the 
number of lattice sites (fc-points) to = 60. 

A. Ground-state properties 

1. Order parameter 

We begin by scanning the parameter space of the ID 
EFKM, in order to detect an EI ground state. The T = 
quantum phase diagram of the ID EFKM has been pre- 
viously explored by the CPMC technique^S, after map- 
ping the EFKM — rewritten in pseudospin variables — into 
a negative U asymmetric Hubbard model with the Zee- 
man term replaced by a chemical potentiaU^. Thereby, 
in terms of our original language, a transition from a 
mixed- valence regime to a nonmixed valence regime was 
observed, which corresponds to the transition from the 
EI to a band insulator. In ID, the EI phase is char- 
acterized by critical excitonic correlations. Surprisingly 
the topology of the ID phase diagram is the same as for 
the 2D and 3D cases, which were studied by CPMGi^ 
(2D) and Hartree-Fock^'^'i^ (2D, 3D) approaches. To 
benchmark the reliability of the PRM, we have compared 
the ID PRM EI band-insulator transition points with 
those obtained by the rather unbiased CPMC method, 
and found excellent agreement. For example, we ob- 
tain = —1-81 ~ ^c.CEMC = —1-80 for ~ —0.3, 
U = 1 (cf. Fig. 1 from Ref/jll. 

Figures [2] and |3] show the onset of the EI phase for 
different values of U {t^ = -0.3 fixed) and {U = 0.8 
fixed), respectively, as the /-electron level is varied (see 
main panels). Obviously, the EI phase emerges above a 
critical Coulomb attraction strength Ud (cf. inset Fig. [2]), 
provided that c and / bands overlap. Moving up the 
e-f level, the top of the / band reaches the bottom of 
the c band at a critical value e^. Then some / elec- 
trons can be transferred into c-band electrons and exciton 
bound-states of /-band holes and c-band electrons may 
form if the Coulomb attraction is sufficiently strong. Re- 
cent Hartree-Fock- and slave-boson-theory-based stud- 
yield a second, upper critical value of the 



les- 



17.18.21 



Coulomb attraction Uc2, such that the EI phase is con- 
fined in between Ud and Uc2- The PRM data, produced 
up to now, will rather not confirm this controversial find- 
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FIG. 2: (Color online) EI order parameter, A, in the ID 
EFKM. The /-electron transfer-integral is fixed to be t-^ — 
—0.3 (all energies are given in units of f^). In what follows, 
we assume e'^ = without loss of generality. The main panel 
gives A as a function of the position of the /-electron level, , 
for different values of U, while the inset shows the variation 
of A with U at £^ — —1.7. The calculations were performed 
at basically zero temperature, T — 10~^. 



ing (see inset Fig. However in the large-C/ limit, the 
numerical calculations are ill-conditioned and tedious, es- 
pecially in ID. So the question whether Uc2 = oo remains 
open. That the appearance of the EI phase is intimately 
connected with the build-up of f-c electron coherence 
and a non-integral /-electron valence is demonstrated by 
the inset of Fig. [31 depicting (n-^) (cf. also the discussion 
of Fig. [5] below) . 



The physical picture developed so far does not change 
if the /-electron band approaches the c-electron band 
from above. That is our results are one and the same 
changing the sign of e-^. We note that in a rather small 
region around e-^ = (symmetric band case), a charge- 
density- wave state is assumed to be the true ground state, 
i.e., the EI phase becomes metastable. Focussing on the 
characterization of the EI phase we will not address this 
issue here. 



2. Band renormalization 

Next we investigate the renormalization of the ID band 
structure. Figure |4] displays the fc-dependence of the 
quasiparticle energies in the EI phase, where open 
and filled symbols correspond to bands having predom- 
inantly c- and /-electron character, respectively. There 
are two features of importance. First, the band structure 
is clearly gapful in the EI phase. Hence, for the half-filled 
band case, the system is insulating. The gap originates 
from c-f electron hybridization, induced by the attrac- 
tive Coulomb interaction. Second, the lower and upper 
quasiparticle bands are narrowed as a result of the elec- 
tronic correlations. While the maximum of the lower (at 
T — Q completely occupied) band is displaced to larger 
/c- values as U increases, the minimum of the upper (at 
T = empty) band moves to smaller fc, accompanied by 
a flattening of the '/-band' near k = 0. At very large U 
(not depicted), the gap is kept (but different in nature), 
because of the extreme Hartree shift, leading to a c-f 
band-splitting. Except of contributions from the band 
narrowing, the Hartree shift to the quasiparticle energies 
is roughly given by the sum of e"^ — and U{{nf) — {rf)) 
(compare Eqs. ([6]). (jiT j) . (j42l) ). 



-0.08- 




FIG. 3: (Color online) EI order parameter in the ID EFKM. 
The main panel gives A as a function of the position of the 
/-electron level for different at U = 0.8. The inset shows 
the corresponding variation of the mean /-electron number. 



3. Momentum distribution function 

Another quantity of interest is the occupation num- 
ber of fermionic states carrying momentum k. For free 
fermions, at T = 0, all states up to the Fermi energy, 
Ep, are occupied, so that the momentum distribution 
function, n(k) = (uk), has a discontinuity at the corre- 
sponding Fermi momentum, kp, where n(k) jumps from 
one to zero . In an interacting Fermi liquid there is still 
a discontinuity, but the jump is less than one. In ID, 
normally Luttinger-liquid behavior emerges, with an es- 
sential power- law singularity at kp. For the insulating 
state, however, n(k) is given by a smooth curve. This 
holds, e.g., for the charge-density-wave ground states of 
IT) t — V and Holstein-type models'^"^, and should also 
be vahd for the EI phase in the ID EFKM. Indeed, the 
momentum distribution functions of c and / electrons, 
rf{k), depicted in Fig. [5] confirm this picture. We see 
that {k) {n'^{k)) monotonously increases (decreases) as 
k varies from fc = to fc = tt. As expected, the drop (up- 
turn) near 'fcp-' softens at larger interactions strengths. 




FIG. 4: (Color online) Renormalized quasiparticle band dis- 
persion of c electrons {E^, open symbols) and / electrons (-E^, 
filled symbols) in the EI phase of the ID EFKM for different 
values of U. Again the 'bare' band structure is parameterized 
by e*^ = 0, = —1.0, — —0.3 (dot-dashed and solid lines). 
Note that the scale of the ordinate is shifted in order to fix 
the Fermi energy at zero energy. 



We have also included in Fig. [5] the variation of dj. = 

iflcj.). This expectation value enters into the equation 
for the order parameter In some sense, it can be 
taken as a measure of the range in /c-space, where c- 
electrons and /-holes are involved in the exciton forma- 
tion and condensation process. Having a {/-driven BCS- 
BEC crossover scenario in the EI phase of the EFKM in 
niind^'^, the broadening of the distribution of d^. with 
increasing U might indicate the condensation of a more 
local two-body EEC-like bound-state out of a BCS-like 
Cooper-pair state. Note that in the BEC-like state the 
Fermi surface plays no role. 



B. Spectral properties 

In this section, we present first results for the pho- 
toemission (PE) spectra in the EI phase of the EFKM. 
The calculated single-particle spectral functions, associ- 
ated with the emission (PE) or injection (inverse PE) of 
an electron with wave vector fc, directly measure the oc- 
cupied and unoccupied densities of single-particle states, 
and therefore are well suited to investigate pairing gaps 
as well. Note that the ARPES spectral functions 



1 /■+°° 



FIG. 5: (Color online) Momentum distribution functions 
(nl) = (4cfc> (black solid hues) and (n{> = (/^/fc) (blue 
dash-dotted lines) for the same model parameters as used in 
Fig. [l] The red dashed lines show the corresponding 'order 
parameter' functions df., see Eq. ((SJ. 



fulfill at T = the frequency sum rule 



(64) 



1 



1 + 6*3" 



A"(k,C^) 



(63) 



where (nj^) is given by Eqs. (gH]) and (cf. Fig.[5|). 

Figure [5] displays the zero-temperature, wave-vector 
and energy resolved single-particle spectral functions, 
A^lkjUj) [see Eqs. (|6T|) . (|62]) ]. for a bare band structure 
parameterized by — —1 (e^ = 0), = —0.3 {f^ = 1). 
For weak Coulomb attraction, U = 0.2 (upper panels) 
we are still in the semi-metallic phase, and consistently 
A''(k,uj) and A-^{k,u!) follow the nearly unrenormalized 
c- and /-band dispersions, respectively. Concomitantly, 
we find a more or less uniform distribution of the spectral 
weight and negligible incoherent contributions (see right- 
hand panels for A'^{k,u!) and Af{k,uj)). When entering 
the EI phase by increasing U to U = 0.6, a gap feature 
develops at the Fermi energy (Fermi momentum), but 
away from that the spectra still show the main charac- 
teristics of the semi-metallic state (cf. both middle pan- 
els). At a still larger value U = 1.2, the gap broadens. 
Most notably, however, is a significant redistribution of 
the spectral weight from the coherent to the incoherent 
part of j4''(fc,a;), with pronounced absorption maxima 
at fc = 0, TT. Of particular importance is a considerable 
admixture of c-electron contributions to the /-electron 
spectrum in an interval between k — and k ^ kp. The 
resulting double peak structure around the Fermi level 
can be considered as an almost /c-independent bound ob- 
ject of c electrons and / holes. 

If we change the location of the /-band by lowering 
the position of the /-electron level, we can achieve that 
the band structure becomes gapful due to the Hartree 
shift even for moderate values of U (assuming A = 




FIG. 6: (Color online) Wave-number resolved photoemission spectra of the ID half-filled EFKM. The c- (left-hand panels) 
and /-electron (right-hand panels) single-particle spectral functions were calculated for several characteristic values of U at 
= —1.0, = —0.3, T = 10~^. The left panels show in each case the total spectra, A^{k,uj), whereas the right panels give 
the 'incoherent' contributions only (second term in Eqs. (|61|) . (|62|) '). The vertical dot-dashed lines mark the chemical potential. 



0). For example, at — —1.7 and — —0.3, a c- 
/ band-splitting (positive Hartree gap A/f, cf. Fig. [J) 
occurs already at Uh ~ 0.94. Such a situation conies 
closer to the Tm[Se,Te] system®. Figure [7] shows the PE 
spectra calculated for these parameter values. Compared 
to the case = —\, ior U = 0.6 the gap is clearly a 
bit more shaped, but the main features of the spectrum 
stay the same. Here we are in the BCS-regime, where 
pairing fluctuations are expected to be small. For U = 
1.2 > Uh: we enter the BEC-regime, where preformed 
pairs acquire quantum coherence (many-body character) 
during the condensation process. Since A ^ the gap 
persists. The distinct incoherent contributions, showing 
up in the yl''(fc,a;) spectra at high energies, are related to 
the dissociation of two-particle bound-states (excitons). 



To visualize more clearly the spectral weight and 
line-shape of the various absorption signals, we pro- 
vide in Fig. [8] a color-map intensity-plot of the single- 
particle spectral functions depicted in Fig.[71 Particularly 
from the [k^uj) data, it appears that the BCS-BEC 
crossover is again marked by a notable admixture of c- 
electron-like contributions leading to new peak-line struc- 
tures in the spectrum (cf. right-hand panels of Fig. [5]). 



V. CONCLUSIONS 

To summarize, we adapted a continuous version of 
the projective renormalization method to the extended 
Falicov-Kimball model to examine the possible existence 
of a collective excitonic insulator phase at the semimetal- 
semiconductor transition. Thereby the PRM approach 
allows us to derive analytical expressions for the EI order 
parameter and various other physical quantities charac- 
terizing the ground-state and spectral properties of the 
model. The self-consistent evaluation of the renormal- 
ization equations yields a stable EI solution for the one- 
dimensional EFKM, at half-filling and zero temperature. 
It therefore confirms previous, straight numerical data 
by constrained path integral Monte Carlo. In particu- 
lar, the phase boundary between the excitonic and band 
insulator agrees even quantitatively with the CPMC re- 
sults. Thus, increasing the Coulomb attraction between 
c-band electrons and /-band holes, the appearance of a 
semimetal-EI transition seems to be settled for the ID 
EFKM. 

Moreover, we present novel results for the single- 
particle c- and /-electron spectral function of the EFKM. 
The calculated photoemission spectra are in evidence of 
a BCS-BEC crossover of the excitonic condensate, trig- 




FIG. 7: (Color online) Single-particle spectral functions for c (left-hand panels) and / electrons (right-hand panels) in ID 
half-filled EFKM with = -1.7, = -0.3, T = 10"^. As in Fig. [51 the total spectra are contrasted to the 'incoherent' 
contributions. 



gered by the Coulomb interaction. Thereby the character 
of the electron-hole pairs changes from the many-body 
bound state associated with the Cooper-type instabil- 
ity (weak-to-intermediate coupling BCS-side) to the two- 
body (tightly-bound exciton) bound-state (intermediate- 
to-strong coupling BEC side) , where Fermi surface effects 
are negligible. Hallmark of the BCS-BEC crossover in 
the quasiparticle spectra is a substantial spectral weight 
transfer from the coherent to the incoherent part of the 
spectrum. A more thorough analysis of the pairing fluc- 
tuations, which are expected to be strongly enhanced 
in the BCS-BEC transition region, e.g. by calculating 



the dynamical pair-susceptibilities within the PRM ap- 
proach, would be a worthwhile goal of forthcoming stud- 
ies. 
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